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Recent evidence for one or more gamma-ray lines at ~ 130 GeV in the Fermi-LAT data from 
the Galactic Center has been interpreted as a hint for dark matter annihilation to Z'y or Hy 
with an annihilation cross section of (cry) ~ 10~ 27 cm 3 s~ 1 . We test this hypothesis by comparing 
synchrotron fluxes due to the electrons and positrons from decay of the Z or the H bosons only 
against radio data from the same region in the Galactic Center. We find that the radio data 
from single-dish telescopes marginally constrain this interpretation of the claimed gamma lines 
for a contracted NFW profile. Already-operational radio telescopes, such as LWA, VLA-Low and 
LOFAR, and future radio telescopes like SKA, are sensitive to annihilation cross sections of the order 
of 10 _28 cm 3 s _1 , and can confirm or rule out this scenario very soon. We discuss the dependencies 
on the dark matter profile, magnetic fields, and background radiation density profiles, and show 
that the constraints are relatively robust for any reasonable assumptions. Independent of the above 
said recent developments, we emphasize that our radio constraints apply to all models where dark 
matter annihilates to Z7 or H-y. 

PACS numbers: 95.35. +d, 98.35.Jk 



I. INTRODUCTION 

The particle identity of dark matter (DM) is one of the 
outstanding puzzles in contemporary physics. In order to 
fully understand the particle properties of dark matter, 
a number of complementary approaches to dark matter 
searches have been adopted. Indirect detection of dark 
matter is a promising technique, in which the products 
of dark matter annihilation are searched for, and gives us 
information about the DM abundance and annihilation 
rate at various astrophysical sites [1-5]. 

Gamma-ray lines from DM self annihilation are 
believed to be a smoking-gun signature, and have been 
investigated in considerable detail [6-11]. Despite the 
relative freedom in DM model-building, if DM self- 
annihilation is to two-body Standard Model final states, 
then gamma-ray line(s) can be produced only via the 
following three channels: (i) XX ~ * 77j (h) XX 
and (iii) \X ~^ Hli where \ denotes the DM, and Z and 
H denote the Z and Higgs boson, respectively. We take 
the mass of the Higgs boson to be 125 GeV [12, 13], and 
allow a heavier DM to annihilate to it. 

Recently, evidence for a gamma-ray line from the 
Galactic Center (GC) has been uncovered in the Fermi- 
LAT data at ~ 130 GeV [14, 15] and this has given 
rise to renewed interest in considering the line signal in 
more detail [16-31]. The Fermi-LAT collaboration has 
not confirmed the signal yet [32], as their most recent 
results have employed a different search strategy, an older 
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data set and background rejection software, and a larger 
search region, making it difficult to compare directly with 
the above claims. 

This statistically significant signal has been tentatively 
interpreted as arising from DM annihilation. Generally 
speaking, the signal requires a DM self annihilation 
cross section of (av) ~ 10 _27 cm 3 s _1 and the Galactic 
DM halo described by a standard NFW, Einasto, or 
a contracted NFW profile. Subsequently, a variety of 
particle physics models have been proposed to explain 
the signal [33-55]. It is also found that the line is off- 
center from the GC by approximately 1.5° [20, 27], which 
requires the center of the DM halo to be displaced from 
the baryonic center. This degree of displacement appears 
reasonable as shown by recent numerical simulations [56] . 

On the other hand, there are arguments against the 
DM origin of the gamma-ray line. There are hints that 
the line is also present in the photons collected from the 
cascades in the Earth's atmosphere, which is a "pure 
background" region [20], although this has been claimed 
to be due to statistical fluctuations [57, 58]. There 
have been claims of the presence of gamma-ray lines at 
the same energy, spatially correlated with some Fermi- 
LAT unassociated sources [24]. However, there are also 
counterclaims that most of these unassociated sources are 
consistent with being standard astrophysical objects such 
as active galactic nuclei or statistical fluctuations [59-61]. 
Furthermore, it remains possible that the GC line signal 
is also of an astrophysical origin [17, 62, 63]. Given the 
arguments in favor of and against the DM origin of this 
signal, this remains a topic of active research. 

If the DM annihilates to two-body Standard Model 
final states, as in (i)-(iii), then we can predict some 
particle physics model-independent consequences. For a 
dominant annihilation (i), i.e., to two photons, there are 



2 



no further interactions of the photons at an appreciable 
level, with all higher-level amplitudes suppressed by at 
least a ps 1/137. However, if the annihilation proceeds 
as in (ii) or (iii), i.e., to a photon and a heavy Standard 
Model boson, the heavy boson decays to other charged 
particles which can have observable consequences. 

The decay of the Z and the H boson produces 
electrons, protons, neutrons, neutrinos, their 
antiparticles, as well as photons as final states. The 
almost featureless spectra of these secondary particles 
poses considerable difficulty in their search above the 
astrophysical backgrounds. Searches for antimatter 
benefit from lower cosmic ray backgrounds, therefore 
one can search for antiprotons and positrons from the 
Z and the H boson. A search for antiprotons from 
these decays constrains several particle physics models 
which can give rise to a gamma-ray line [21], whereas 
the preexisting unaccounted excess in positrons [64, 65] 
makes a positron search ambiguous. Neutrinos could, 
in principle, be used to distinguish between all three 
final states, but achieved or projected sensitivities in the 
range (crv) ~ (1CT 22 cm 3 s" 1 - 1CT 23 cm 3 s" 1 ) [66-69] 
will not be able to probe the claimed signal. Secondary 
photons that are produced in the decay of the Z or 
the H boson, or in other DM annihilation channels, 
also constrain these scenarios [21-23], and there are 
ongoing efforts to confirm this 130 GeV line with future 
detectors [25, 30]. 

In this paper, we ask the questions - If the 130 GeV 
signal is indeed from DM annihilation to Z-f or Hj, what 
other consequences are guaranteed? Can we use these 
consequences to test this signal? Synchrotron radiation 
from products of DM annihilation has been argued to 
provide strong constraints for many DM annihilation 
channels and scenarios [70-86]. Thus, following these 
promising leads, we explore our question by calculating 
the synchrotron radiation in the GC due to the electrons 
and positrons from Z or H decays only, and comparing 
it to existing data from radio telescopes. 

We first take a very conservative approach, where 
we compare the DM-induced synchrotron fluxes to the 
total measured radio flux at 330 MHz in a relatively 
large region around the GC, and determine that DM 
annihilation cross sections to these channels cannot be 
more than (av) ~ 10 _25 cm 3 s . However, this approach 
is overly conservative, as the synchrotron fluxes in the 
GC are modeled accurately with known astrophysics. 
We argue that the flux due to dark matter must not 
exceed the uncertainties on the modeled radio fluxes, 
which provides us with a constraint (av) ~ 10 cm s . 
Constraints obtained by comparing fluxes predicted in 
smaller regions of interest and upper limits at 408 MHz 
imply (av) ~ 10~ 27 cm 3 s -1 , and are already in mild 
tension with the 130 GeV line. We forecast that the 
sensitivity can be improved to ~ 10 _28 cm 3 s _1 with a few 
hours of observation of the GC at 80 MHz with LWA, the 
Long Wavelength Array, and at 200 MHz with LOFAR, 
the LOw-Frequency ARray for radio astronomy, allowing 



us to constrain interpretations of the 130 GeV line signal 
in the Fermi-LAT data in terms of DM annihilation to 
Zj or Hj. Although this is the main motivation for our 
present work, the radio constraints we derive are valid 
regardless of whether this claimed 130 GeV line signal 
survives further scrutiny or not. These constraints will 
continue to apply to any future interpretations of gamma- 
ray lines at the GC in terms of DM annihilation. 

Note that, these sensitivities readily probe the cross 
section that explains the tentative 130 GeV line signal. 
More generally, we expect these sensitivities to be 
able to probe many of the models, not necessarily 
supersymmetric, that could explain this signal. We also 
emphasize that since we are looking for the synchrotron 
radiation from the electrons and positrons produced in 
the decays of the Z or the H boson only, our constraints 
are independent of the underlying DM particle physics 
model. In these two ways, our work is complementary 
to Ref. [22]. The results here are of course affected 
by astrophysical uncertainties, e.g., dark matter density 
profile, magnetic fields, interstellar radiation energy 
density, and proton density in the Galaxy, and by taking 
a range of different values for them we try to understand 
their impact. 

The rest of the paper is organized as follows. In Sec. II 
we discuss the radio data that we use for obtaining our 
constraints, and the theoretical framework for calculating 
the flux densities from synchrotron radiation by DM 
annihilation products. In Sec. Ill we furnish and justify 
the astrophysical inputs, i.e., DM density, magnetic 
fields, and radiation density in the Galaxy, that we use 
for our calculations. In Sec. IV we show the predicted flux 
densities for benchmark DM annihilation cross sections, 
and provide constraints on the DM annihilation cross 
section as a function of DM mass for the channels 
XX Z"{ and XX ~> H-f, and conclude in Sec. V. 



II. EXPERIMENTAL DATA AND 
THEORETICAL FRAMEWORK 

A. Radio data and regions of interest 

We use radio data at two frequencies, 330 MHz and 
408 MHz, to obtain the limits on DM self annihilation 
cross section from near the GC. We use the projected 
radio sensitivity of the LWA telescope at 80 MHz [87] 
to predict the future sensitivity on DM self annihilation 
cross section to gamma-ray lines that can be probed by 
radio data. We calculate the synchrotron flux density at 
a region offset from the GC in the 200 MHz band, which 
is an operating frequency band for the LOFAR telescope. 

We consider both the LWA and LOFAR telescopes for 
two reasons. First, the geographical location of LOFAR is 
not ideal to observe the GC, but it can measure the radio 
flux away from the GC to derive useful constraints on 
dark matter properties. Second, LWA is in a much better 
location to study the GC, but the GC can be opaque at 
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the frequencies the LWA will operate in. The redeeming 
factor is that the absorptive nature of the GC at these 
frequencies is patchy and there exists regions that are 
transparent [88-91]. This argues for region optimization 
depending on the observed patchiness and sizes. 

1. ROI-2° : Region of interest 2° around GC 

The radio measurements in the 330 MHz band by 
the Green Bank telescope are available in a 6° x 2° 
region around the GC [92], and provides us our first 
region of interest (ROI-2°). We approximate this region 
to be a circle with a radius 0.034 radians 2°), for 
simplicity. Thus, we approximately match the area of 
the region of observation in Ref . [92] , but the shape is 
different. We assume that this difference will not change 
our results significantly. In Ref. [92], the authors present 
an astrophysical model to explain the data, so we use 
the uncertainty in their measurement at 330 MHz, i.e., 
0.05 x 18000 Jy = 900 Jy, to obtain our limits on the 
self annihilation cross section as a function of the DM 
mass. The authors in Ref. [92] also use radio data at 
higher frequencies to construct a GC model. Comparing 
our calculated synchrotron flux density with the errors in 
their measurement at every measured frequency we find 
that the most constraining limit on DM self annihilation 
cross section comes from the lowest frequency band 
(330 MHz) and hence we only use the uncertainty in the 
measurement at 330 MHz to constrain DM properties. 

We will also use this ROI to obtain our projected 
sensitivity on DM particle properties using the future 
measurement around the GC by the LWA telescope. 

2. ROI-4" ■' Region of interest 4" around GC 

The radio measurement in the 408 MHz band by the 
Jodrell Bank telescope [93] in a region 4" around the GC 
provides us with our second region of interest (ROI-4"). 
At this frequency, the region of interest is a circle of 
angular radius 4" and the upper limit on the radio flux is 
50mJy. This region is significantly smaller than ROI-2°, 
and as we will show, is affected differently to our input 
parameters. Thus it provides a complementary site to 
testing DM properties. 

3. ROI-away: Region of interest away from the GC 

We also calculate the synchrotron flux within an 
angular cone of radius 1°, at angles 1.5° and 10° away 
from the GC. We calculate the synchrotron flux at regions 
away from the GC at 80 MHz, which is an operating 
frequency band for the LWA telescope, and at 200 MHz, 
which is an operating frequency band for the LOFAR 
telescope. We calculate how the synchrotron flux varies 
with mass of the DM, for a given (av). The advantage 



of measuring the synchrotron flux away from the GC is 
that the synchrotron flux depends less on the assumed 
DM profile and has much smaller backgrounds. Ideally 
the best sensitivity to DM properties will be found if the 
radio measurement is done in a radio "cold spot" , where 
no known radio sources exist. On the other hand, due to 
smaller DM density, the synchrotron flux falls away from 
the GC. This disadvantage is partially mitigated by the 
excellent sensitivity of the present and upcoming radio 
telescopes like LWA, LOFAR and SKA. 

Finally, anticipating future radio measurements near 
the GC, we also estimate the constraints that can be 
obtained on DM self annihilation channels that produce 
a gamma-ray line using the projected sensitivity of LWA. 
We very conservatively assume that LWA can reach a 
background subtracted sensitivity of 10 Jy at 80 MHz 
in a circular region of radius 2° around the GC. After 
the measurement of the radio flux near the GC one 
has to model the synchrotron flux due to expected 
astrophysical processes and then use the uncertainty 
in that measurement to constrain the DM particle 
properties (as we have done for the 330 MHz band). 
Although we will only present our calculated synchrotron 
fluxes in regions away from the GC, these can also be used 
to measure the DM properties. We remind the reader 
that although the GC is generally opaque to frequencies 
< 100 MHz, the absorption is patchy and the patchiness 
can be used to do the GC radio measurements [88-91]. 



B. Theoretical framework 

To calculate the synchrotron flux from DM self 
annihilation products, in principle, we need to solve 
the time- independent diffusion equation for the produced 
electrons and positrons [73] 

K(E)V 2 n e (E,r) + A [b(E,r)n e (E,r)} 

= -S(E,t), (1) 

where n e (E, r) is the electron density spectrum per unit 
energy interval, K(E) is the diffusion coefficient, b(E, r) 
is the energy-loss rate and S(E, r) is the source injection 
spectrum of the electrons. Here E denotes the energy of 
the electron and r denotes the position of the electron. 

For ROI-2°, it can be shown that we are in a regime 
where the GeV electrons will travel only ~ 30 pc [75, 94] 
during their cooling lifetime. Since this length is much 
smaller than the length associated with ROI-2°, we 
conclude that diffusion will have a small impact on 
our results. For ROI-4", due to the presence of very 
high magnetic fields near the GC (see Sec. IIIB), and 
consequently high energy loss rates, the diffusion length 
(1(E)) ~ ^KE/b, is very small and diffusion can be 
safely neglected. However, for increased precision, one 
may in future improve our results by performing a more 
detailed numerical study along the lines of [86] . 



4 



These electrons and positrons then lose energy via the 
synchrotron process, the inverse Compton process, and 
the bremsstrahlung process. For our purposes, the total 
energy loss rate is given by 

b(E,r) = b syuc (E,r) + b IC (E,r) + (25, r). (2) 

Ionization energy loss is important for electrons and 
positrons only at lower energies than are considered in 
this work. 

Synchrotron energy losses are due to the interaction of 
the electron and the positron with the Galactic magnetic 
field. The energy loss rate due to synchrotron process is 
given by (all formulae in this section are in SI units) [73] 
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where ot = e 4 / (6irelmlc 4 ) is the Thompson scattering 
cross section [95], J7 mag is the magnetic energy density, 
7 = E/m e is the Lorentz factor, and /3 
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The photons emitted because of synchrotron energy loss 
is generally in the radio band. 

Inverse Compton energy losses are caused by the 
up-scattering of the photons in the GC region (which 
is mainly composed of the CMB and the background 
starlight) by the more energetic electrons and the 
positrons. The energy loss rate due to the inverse 
Compton process is given by [73] 
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where C/ ra d is the radiation energy density. The photons 
from the CMB and the background starlight is generally 
up scattered to gamma-ray energies by the energetic 
electrons and positrons from the decays of the Z-boson 
and the Higgs boson. 

Bremsstrahlung losses are caused by the emission of 
gamma-ray photons by the electrons and positrons due 
to their interaction with the nuclei in the Galaxy. The 
energy loss rate for this process is given by the Bethe- 
Heitler formula [96]. We assume that the hydrogen are 
the dominant nuclei present in the Galaxy. The energy 
loss rate due to bremsstrahlung is given by 



dE 
~dt 
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where we use nn ~ 3 cm 3 as the number density of 
hydrogen nuclei in the interstellar matter in the Galaxy. 

The source term is due to the particle injection from 



DM self annihilation 
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where m x denotes the mass of the DM particle, p x { r ) 
denotes the DM density distribution. dN e /dE denotes 
the number of electrons and positrons from the decay of 
the Z or the H boson per unit energy interval, which we 
calculate using PYTHIA [97]. 

Collecting the above mentioned inputs and in the 
no diffusion limit [75], we can write electron density 
spectrum per unit energy interval as 



n e (E,r) = - 



J™ x dE' S(E',r) 
b(E,r) 



(7) 



1. Flux density in ROI-2° 

The synchrotron power density per unit frequency from 
a spectrum of electrons and positrons is given by [96, 98] 



L„(r)= J dE n e {E,r) 

1 v / 3e 3 S(r) / v 



< - A — I — / dx K 5/3 (x) \\ , (8) 
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where the critical frequency, v ci is given by 

3eE 2 B ( E \ 2 / B\ 

^=4^ = 16MHz (Gev) U)' (9) 



and K 5 / 3 (x) denotes the modified Bessel function of 
order 5/3. 

The synchrotron radiation flux density is given by 



(10) 



where I denotes the line of sight distance and fi is the 
angular area of the region of interest. 

We have verified that the synchrotron self-absorption 
is unimportant for these parameters which is generally a 
problem near the GC at frequencies below approximately 
100 MHz but the absorption regions are patchy [88-91]. 



2. Flux density in ROI-4" 

For the smaller region of interest of radius 4" around 
the GC, i.e., ROI-4", we follow the method presented in 
Ref. [72], which is dependent on the morphology of the 
magnetic field near the GC black hole. Due to the strong 
magnetic fields in this region, we assume that the energy 
loss of the electron is dominated by the synchrotron 
energy losses. In this case, we approximate 
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The synchrotron flux density in this case is given 
by [72] 



F u = 
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where the first integral is over the volume of observation 
and the second integral counts the number of particles 
above a certain energy E. The value of E in this case 
can be found by using Eqs. (9) and (11) and is given by 



E = 433 MeV, 



Y v _ \ ( 
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3. Flux density in ROI-away 

We calculate the synchrotron flux in ROI-away in the 
same way as we do in ROI-2°. To account for the fact 
that we are now calculating the synchrotron flux away 
from the GC, we do make some modifications to our input 
of the interstellar proton density and the interstellar 
radiation density. For simplicity, we take the number 
density of hydrogen nuclei in the interstellar medium of 
our Galaxy to be 3 cm -3 . We take the variation of the 
radiation energy density following Ref. [73] . More details 
about the variation of the radiation energy density is 
given in Sec. Ill C. 



III. ASTROPHYSICAL INPUTS FOR 
CALCULATIONS 

A. DM density profiles 

One of the major unknowns near the GC is the DM 
density profile. Almost all simulations agree on the radial 
dependence of the DM density profile at large radii from 
the GC, p(r) ~ r~ 3 . However, due to limited numerical 
resolution and the complicated astrophysics at the GC, 
the simulations disagree on the shape of the density 
profile at small radii. 

Observations of elliptical galaxies [99], early- type 
galaxies [100], M31 [101], and M84 [102] prefer a cuspy 
profile [103] in contrast to dwarf galaxies which prefer 
a cored profile [104]. Given that M31 is a Milky Way 
like galaxy, we assume that the DM density profile in the 
Milky Way is not cored and hence we do not consider 
the cored isothermal profile in our work. In general, 
the constraints from the indirect detection searches are 
especially weak for a cored isothermal profile [72, 75]. 

In this work, we use three different DM profiles 
which provide reasonable constraints from the radio 
measurements at the GC. The dark matter density at 
the solar radius has a value of 0.3 ± 0.1 GeV/cm 3 [105]. 
For concreteness, in this work, we take the benchmark 
value to be 0.4 GeV/cm 3 [106, 107]. Note that local 
dark matter density taken in the papers which discuss 



the presence of the 130 GeV line at the GC is also 
0.4 GeV/cm 3 (see, e.g., Refs. [14, 16]). 

The Einasto dark matter profile [108-111], 



pEin(r) 



0.08 GeV cm" 



exp 




20 kpc 



0.17 



- 1 



(14) 



is the least cuspy of all the DM profiles considered in 
this work, and hence we expect this profile to produce 
least amount of synchrotron radiation, especially when 
we consider the synchrotron radiation from a region very 
near the GC. 

We then consider the standard NFW profile [112] 
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The cuspy nature of this DM density profile will ensure 
that we get a larger synchrotron radiation flux than what 
we expect from the Einasto profile when we consider 
observation from a region very near to the GC. 
We finally consider a contracted NFW profile 



Peon. NFW ( r ) 



0.29 GeV cm 
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20 kpc 



1.15 
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The steeper inner slope in this case can be due to either a 
GC black hole [113], or due to adiabatic contraction due 
to the presence of baryonic matter at the GC [114-117], 
which have been supported by more recent numerical 
simulations [103]. 

These DM profiles are shown in the left panel in Fig. 1. 
It is evident from the figure that, at small radii, the 
contracted NFW profile has the steepest slope, and the 
Einasto profile has the shallowest slope of all the three 
DM profiles considered in this work. From the figure, one 
can also infer that the DM density profiles have almost 
the same shape at large distances from the GC. 



B. Magnetic fields 

The GC magnetic field has both a regular and a 
turbulent component. For both the components, the 
normalization and the radial profile is not understood 
very well. In particular, the magnetic field amplitude 
near the GC is uncertain, with measured estimates 
spanning a range of some two orders of magnitude 
between 10 pG [118] and 10 3 pG on scales of a few 
~ 100 pc [119]. In order to account for the uncertainty 
in the magnetic field, we adopt several configurations. 
In all cases, we initially fix the normalization to a value 
Bq — 6 pG at the solar system radius (r = 8.5 kpc). 
This is mid-range among the various estimates of Bq 
which span between 3 pG and 10 pG [120-122]. We will 
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FIG. 1. Galactic dark matter density (left panel) and magnetic field profiles (middle and right panels). (Left) We show the 
various DM profiles used in this work: the Einasto profile (14), the NFW profile (15) and the contracted NFW profile (16). 
(Middle) We show the magnetic field used for calculating the synchrotron flux in ROI-2°, i.e., a region of angular radius 2° 
around the GC. The constant magnetic field has a value of 6/iG everywhere in the Galaxy. The exponential magnetic field is 
given in Eq. (17) and has a value of 6/iG at the solar radius. (Right) We show the magnetic field used for calculating the 
synchrotron flux in ROI-4", i.e., in a region of angular radius 4" around the GC. The "equipartition" magnetic field is given 
by Eq. (18) and the "cored" magnetic field is given by Eq. (20). Both the fields have been normalized to have a value of 6/iG 
at the solar radius. 



discuss how our results scale with the different values of To estimate the uncertainty due to the shape of the 
Bq when we present our results. magnetic field profile we also adopt the extreme case of 

a constant magnetic field of value Bq everywhere [78]. 

Both of these are shown on the middle panel of Figure 1 . 

1. Magnetic fields in RO I -2° 



For ROI-2°, i.e., a circular region with radius 2° around 
the GC (distance scale ~ 200 pc for r Q = 8.5 kpc), we 
consider two different magnetic field radial profiles. The 
first is the spherically symmetric exponential profile, 

B(r)=B (a e*p(-^^\ , (17) 

where r is the distance from the GC, and R m = 14 kpc is 
the scale radius. Our choice of R m follows from Ref. [74], 
where we adopt their Galactic magnetic field model 
"GMF I" and their best-fit propagation parameters. We 
add that we neglect the z-dependence of the magnetic 
field which is only weakly constrained by data and remain 
highly uncertain. Using this conservative magnetic field, 
the magnetic field at a radius of 2° is ps 11 /iG. Although 
this value is within the range of estimates of the magnetic 
field in the GC, it is closer to the lower range. In 
addition, it does not obey the lower limit of 50 /iG on 
scales of 400 pc presented in Ref. [92]. However, given 
the uncertainties in the astrophysical and propagation 
quantities, we do not consider this difference significant. 
For example, if we adopt instead Bq = 10 /iG and the 
"MAX" propagation parameters of Ref. [74] , we obtain 
R m fa 8.5 kpc and a magnetic field at 2° of pa 27 /iG. To 
estimate the impact of the normalization of the magnetic 
field, we will also show our results for two extreme values 
for Bq, i.e., 3/iG and 10 /iG. 



2. Magnetic fields in ROI-4" 

For ROI-4", i.e., a region with radius 4" around the 
GC (distance scale ~ 0.2 pc for r Q = 8.5 kpc), we need to 
take into account the influence of the supermassive black 
hole (SMBH) at the GC. The presence of the SMBH sets 
two length-scales: the Schwarzschild radius, i?BH ~ 1-2 x 
10 12 (A//4.3 x 10 6 M Q )cm, and the radius i? acc « 0.04 pc 
within which the free-fall velocity due to the gravity of 
the SMBH v — — cy Rbh/ 7 " is larger than the random 
Galactic motion, ~ 10~ 3 c. In other words, the region 
r < i? acc defines the accretion region. 

We adopt the "equipartition model" for the Galactic 
magnetic field, described by various authors [e.g., 72]. 
In this model, the SMBH accretes matter from a radius 
of -R acc , and the magnetic field in the accretion flow 
achieves its equipartition with the kinetic pressure, i.e., 
.B 2 (r)/(2/io) = p(r)v 2 (r)/2. For a constant mass 
accretion rate, M, one obtains p(r) — M / Airr 2 v{r) cx 
r~ 3 / 2 , and thus B(r) cx r~ 5 / 4 . Outside of i? acc , the 
conservation of magnetic flux is assumed, yielding B(r) cx 
r~ 2 . Thus, the equipartition magnetic field is given by 

( B acc {r/R acc )~ 5 / 4 r < R acc 
B cq (r) = I B acc {r/R acc )- 2 i? acc < r < i? flux (18) 
[ Bq i?fl ux < r , 
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where i? acc « 0.04 pc, R Rux - 100i? acc , and 

B »= 7 - 9m K«^k;) 1,4 (i«j^)' (19) 

for typical values of B & . 

We also consider a variant of the equipartition model, 
where the inner magnetic field is kept smaller because 
equipartition is prevented somehow. This may occur 
if, for example, magnetic field dissipation occurs by 
reconnection in the turbulent accretion flow (see, e.g., 
[123] and references therein). Since the details of 
dissipation are not well understood, we conservatively 
adopt a constant magnetic field throughout the accretion 
region, namely, 

{B a cc ^ — Race 

B acc {r/R acc )- 2 R acc <r< i? flux (20) 
B Q i? flux < r . 

We call this the cored magnetic field. These are shown 
on the right panel of Figure 1. 



3. Magnetic fields in ROI-away 

While calculating the magnetic field in a region offset 
from the Galactic Center, we assume the exponential 
magnetic field structure as in Eq. (17). 



where we denote the CMB energy density as Uqmb and 
the stellar radiation energy density by U s t e iiar- The 
vertical distance from the plane of the Galaxy is denoted 
by z and the radial distance from the center of the Galaxy 
is denoted by r. 

IV. RESULTS 

In this section, we first present our results for the 
synchrotron fluxes from DM annihilation products, and 
discuss the expected systematic uncertainty due to 
incomplete knowledge of the DM and magnetic field 
profiles. In all the plots, we have taken (av) — 
10 -26 cm 3 s _1 unless otherwise mentioned. We only 
assume the DM self annihilation channels \X ~~ * or 
XX ~~ Hi t° present a completely particle physics model 
independent result. We then compare the expected flux 
with available/projected radio data from the GC to arrive 
at constraints/sensitivities on the DM annihilation cross 
section. We perform this exercise for three frequency 
bands (330 MHz, 408 MHz, and 80 MHz) in two different 
regions of interest around the GC. We also calculate the 
synchrotron flux due to DM annihilation using the above 
mentioned parameters in a region of radius 1° at an angle 
10° away from the GC. Although we do not use the 
region offset to the GC to derive any constraints on DM 
properties, we remind the reader that strong constraints 
can also be obtained from radio observation at a region 
away from the GC. 



C. Radiation energy density 

The radiation energy density is the sum of the energy 
density of the CMB photons and the energy density of 
the background starlight photons. The energy density of 
the CMB photons is 0.3 eV cm" 3 . The radiation energy 
density due to the background starlight varies with 
position in the Galaxy from 1 eV cm" 3 to 10 eV cm" 3 . 
Conservatively, when we calculate the synchrotron flux 
in a region near the GC, i.e., in ROI-2° and ROI-4" , 
we take the background starlight energy density to have 
a constant value of 9 eV cm" 3 . Hence we use the total 
radiation field energy density as t/ rac i = 9 eV cm" 3 in this 
work while calculating the synchrotron flux in a region 
near the GC. We also check our results by taking a much 
smaller radiation energy density of 0.9 eV/cm 3 and find 
that using this lower value of the radiation field energy 
density improves our constraints by a factor of ~ 2-3. 

When we calculate the synchrotron flux in a region 
away from the GC, ROI-away, we follow the radiation 
energy density parametrization in Ref. [73] which uses the 
results in Ref. [124]. The energy density at any given 
position in the Galaxy is [73] 

77 t \ ^stellar (4 kpc,z) 

E/rad(f, Z) = — — , E/ St ellar(r, 0) + U CMB , ( 21 ) 

(^stellar (4 KpC, UJ 



A. Synchrotron flux density at the GC 

1. Flux density in ROI-S 

We calculated the synchrotron flux according to the 
prescription and inputs presented in Sec. Ill and plot, 
in the top and bottom panels of Fig. 2, the flux density 
due to synchrotron radiation from DM self annihilation 
products against mass of the DM for a region of 
radius 2° around the GC in the 330 MHz and 80 MHz 
radio band respectively. The magnetic field used is a 
spherically symmetric exponential magnetic field, taken 
from Eq. (17) with the magnetic field at the solar 
radius normalized to be 6 p,G. The uncertainty in the 
measurement of the synchrotron flux density at 330 MHz 
band in this region around the GC at this frequency 
band, 900 Jy [92], is also shown in the plot and is used 
to obtain our constraints on the DM particle properties. 

For both the frequency bands, the synchrotron flux 
from DM annihilation products is maximum for the 
contracted NFW profile. This is expected because the 
signal is proportional to p 2 and a larger p increases the 
signal at the GC. When the region of interest is fairly 
large, e.g., ROI-2°, the Einasto profile is predicted to 
lead to more annihilation than the standard NFW profile. 
At such a large distance from the GC, the synchrotron 
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FIG. 2. Prediction of the synchrotron flux density 2° around the GC, against mass of the DM. In all the plots, we set as 
a benchmark DM annihilation cross section {av} = 10 -26 cm 3 s _1 , and consider three different DM profiles: Einasto profile 
in Eq. (14), NFW profile in Eq. (15) and the contracted NFW profile in Eq. (16). (Left) X X ~> Zj. (Right) XX -> #7- (Top) 
Results for 330 MHz. We also show the uncertainty in the measurement which is used to derive our constraints in this frequency 
band. (Bottom) Results for 80 MHz. For both the frequency bands, we use the exponential magnetic field in Eq. (17). 



energy density only varies by a factor of a few for different 
DM profiles, demonstrating the relative robustness of 
these results. 



2. Flux density in ROI-4" 

We calculated the flux densities and plot it against the 
mass of the DM in the top and bottom rows in Fig. 3 for 
the 408 MHz and 200 MHz radio band respectively. The 
changes are that the region of interest is now a circular 
region of radius 4" around the GC and the frequency of 
radio observations is taken to be 408 MHz. The magnetic 
field used is the equipartition magnetic field, taken from 



Eq. (18). We get very similar results (differences of less 
than 1 mjy) if we use the cored magnetic field, as given 
in Eq. (20). A different value of the magnetic field at 
the solar radius (within the range 3 /j,G to 10 /iG) does 
not change the value of the synchrotron flux density by 
more than a factor of two. The upper limit on the 
synchrotron flux density in this region around the GC 
at this frequency, 50 mjy [93], is also shown in the top 
plot and is used to obtain our constraints on the DM 
particle properties. 

As expected, the synchrotron flux from DM 
annihilation products is maximum for the contracted 
NFW profile. However, in contrast to the above we 
find that for smaller regions of interest, e.g., ROI-4" the 
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FIG. 3. Prediction of the synchrotron flux density 4" around the GC, against mass of the DM. The DM annihilation cross 
section and the profiles are the same as in Fig. 2. (Left) \X ~^ (R-ight) XX ~~ > Hy. (Top) Results for 408 MHz. (Bottom) 
Results for 200 MHz. We use the equipartition magnetic fields in Eq. (18) for both these frequency bands. 



cuspiness of NFW profile at lower radii leads to larger 
fluxes than from the Einasto profile. Note however, that 
for such small regions of interest around the GC, the 
flux varies by orders of magnitude for the different DM 
profiles. 



3. Flux density in ROI-away 

We calculated the synchrotron flux according to the 
prescription and inputs given in Sec. Ill in ROI-away 
and plot some representative results in Fig. 4. The 
top and bottom panels show the synchrotron flux in a 
circular region of radius 1° at 10° away from the GC 
for the 200 MHz and 80 MHz radio band respectively. 
The magnetic field used is the exponential magnetic field, 



taken from Eq. (17) with the magnetic field at the solar 
radius normalized to be 6 fiG. We also take into account 
the variation of the radiation density with distance from 
the GC while calculating the synchrotron fluxes following 
the parametrization in Eq. (21). For a given angle 9 away 
from the GC, we calculate the value of z and then use the 
radiation energy density [/stellar (0, z) in our calculations. 
Although this is a conservative approximation, we expect 
that a full calculation will given similar results. 

The disadvantage of the synchrotron flux decreasing as 
one makes a measurement away from the GC is mitigated 
by the fact that the flux depends less strongly on the 
assumed DM profile. Due to the excellent sensitivity of 
present generation radio telescopes like LWA and LOFAR 
and even better sensitivity of near future radio telescopes 
like SKA very robust limits on DM properties can be 
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FIG. 4. Prediction of the synchrotron flux density at a region if radius 1° at 10° away from the GC, against mass of the DM. 
The DM annihilation cross section and the profiles are the same as in Fig. 2. (Top) Results for the 200 MHz radio band. 
(Bottom) Results for the 80 MHz radio band. We use the exponential magnetic fields in Eq. (17). 



obtained from radio measurements away from the GC. 
In particular, as mentioned earlier, if the measurement is 
done in a radio cold spot then modeling the astrophysical 
backgrounds will also be easier to find the putative radio 
signal of DM annihilation. 



4- Variation of the synchrotron flux with angle away from 
the GC 

We plot our calculated synchrotron flux density against 
angle away in the GC for the 200 MHz band and for the 
XX channel in Fig. 5. We assume the exponential 

magnetic field as given in Eq. (17). We calculate our 
synchrotron fluxes in a region of radius 1° at the specified 
angle away from the GC. We take into account the 



variation in the radiation field energy density following 
the prescription in Sec. IV A 3. As can be seen from 
the plots, the synchrotron flux falls off by an order of 
magnitude as the angle away from the GC increases from 
10° to 50°. A very similar variation of the synchrotron 
flux away from the GC is obtained for the XX ^ Hj 
channel and in the 80 MHz radio band. 



B. Sensitivity to magnetic fields 

Now, we explore the sensitivity of the predicted 
synchrotron fluxes to the normalization and shape of the 
Galactic magnetic field profile. We remind the reader 
that we have used 6 fj,G as our benchmark value of the 
magnetic field at the solar radius for all our calculations 
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FIG. 5. Prediction of the synchrotron flux density against 
angle away from the GC for a DM mass of 100 GeV at 
200 MHz. The region of interest is a circular region of radius 
1° at the specified angle away from the GC. We take the 
NFW DM profile, the same annihilation cross-section, and 
magnetic field as in Fig. 4. The variation of the synchrotron 
flux density with angle is very similar for all the DM profiles 
and annihilation channels considered in this work. 



FIG. 6. Prediction of the synchrotron flux density vs. mass 
of the DM at 330 MHz, in a region 2° around the GC. The 
magnetic field used is exponential, as in Eq. (17), and a 
constant magnetic field. We use two different values of the 
magnetic fields at the solar radius: 3 pG and 10 fj,G. We use 
the NFW profile and the same annihilation cross section as 
in Fig. 2. The variation is similar for all the DM profiles and 
annihilation channels considered in the text. 



presented in the other sections. The DM is assumed to 
have a standard NFW profile (Eq. (15)) for the plot in 
this section. 

To understand the impact of the normalization of the 
Galactic magnetic field on the synchrotron flux density 
due to dark matter annihilation, in Fig. 6 we plot the 
synchrotron flux due to two different values of the 
Galactic magnetic field at the solar radius: 3 /jG and 
10 fiG. It is seen that varying the normalization of the 
Galactic magnetic field can change the synchrotron flux 
density by a factor of a few for both the exponential 
magnetic field profile and the constant magnetic field 
profile. 

We show the impact of two different magnetic field 
profiles: the exponential profile, Eq. (17), and the 
constant magnetic field profile for the 330 MHz band. 
For a given normalization of the magnetic field profile 
at the solar radius, we see that the flux due to the 
exponential magnetic field profile is always larger than 
the flux due to the constant magnetic field profile for the 
DM annihilation channel, \X ~^ The result is similar 
for the annihilation channel \X ~ * Hj. 

The overall uncertainty in the normalization and the 
shape of the Galactic magnetic field can lead to a 
difference of at most an order of magnitude in the 
predicted synchrotron flux in the 2° around the GC 
at 330 MHz. We have checked that the variation in 
the synchrotron flux density with the normalization and 
shape of the Galactic magnetic field is similar for the 
other DM profiles considered in this work. 



For the region of angular radius 4" around the Galactic 
Center, the difference in the synchrotron flux density 
is less than a factor of two for both the equipartition, 
Eq. (18) and the cored magnetic field profile, Eq. (20), 
for a given DM profile. Again, the uncertainty due to 
incomplete knowledge of magnetic fields can lead to at 
most an order of magnitude changes in the predicted 
synchrotron fluxes. 



C. Constraints on {av}-m 

1. Constraint from the measurement at 330 MHz 

In the region with radius 2° around the GC, data in 
the 330 MHz radio band is presented in Ref . [92] . We 
compare our prediction of the synchrotron flux from 
products of DM annihilation, and demand that the 
radio is not over-saturated by the DM-induced fluxes. 
This gives us a constraint on the DM annihilation 
cross section (<jv) ~ 10 _25 cm 3 s _1 . This is overly 
conservative, as there are known astrophysical sources 
that produce most of the observed synchrotron radiation. 
The astrophysical model presented in Ref. [92] suggests 
that with present level of uncertainty, at most 5% of the 
flux (< 900 Jy) could come from unknown sources. This 
gives a much stronger constraint (cry) ~ 10~ 26 cm 3 s -1 . 
This constraint on the (av)-m plane that can be derived 
from the radio flux measurement at 330 MHz for a 
circular region of radius 2° is plotted in the top panel of 
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FIG. 7. Constraints obtained in the (av) vs. m plane for different annihilation channels, different frequency bands, and different 
regions of observations. The DM profiles are the same as in Fig. 2. (Left) \X ~^ (Right) \X ~^ Hi- (Top) Results for 
330 MHz using the exponential magnetic field in Eq. (17). (Middle) Results for 408 MHz using the equipartition magnetic field 
in Eq. (18). The DM mass and (av) preferred for the 130 GeV gamma-ray line is shown by the thin box (see text for more 
details). (Bottom) The sensitivities at 80 MHz using the exponential magnetic field in Eq. (17). The shaded region shows 
approximately the (av) favored by the recent 130 GeV gamma-ray line signal for the NFW profile (green), Einasto profile (red) 
and contracted NFW profile (blue) [14]. The favored region in masses is same as the plot above, and are staggered in this 
bottom plot for better legibility. We also show the updated total thermal relic cross section for reference [125]. 
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Fig. 7. We only show the constraints that can be obtained 
in this radio band by using the exponential magnetic field 
given in Eq. (17), with a normalization of 6/iG. 

For both the DM self annihilation channels XX ~~ * %1 
and XX ~^ Hi: we see that the contracted NFW profile 
gives the most constraining limit. Since the gamma- 
ray line prefers a cross section (av) ~ 10 _27 cm 3 s _1 for 
all the three profiles [14], it can be concluded that the 
existing data at this frequency is not able to constrain the 
line signal independent of a DM particle physics model. 
However, since the present constraints are only an order 
of magnitude away from the DM self annihilation cross 
section preferred by the 130 GeV signal, a future radio 
measurement near the GC can be used to either constrain 
or confirm its presence at the GC. 



2. Constraint from the measurement at 408 MHz 

The upper-limit on the synchrotron flux at 408 MHz 
found by Ref. [93] allows us to impose much stronger 
constraints than above. The procedure that we follow 
is similar to above - we compare the predicted fluxes 
with the existing upper limit, and demand that the 
DM annihilation not produce a flux larger than what is 
already constrained. This constraint in the (av)-m plane 
that can be derived from the radio flux measurement at 
408 MHz for a circular region of radius 4" is plotted in 
the middle row of Fig. 7. 

We also show the dark matter mass and self- 
annihilation cross section preferred for the 130 GeV 
gamma-ray line by the thin shaded box. For the 
annihilation to Zj/ Hj, the gamma-ray energy is 

given by E 1 = m x (l — m z ^ h / Arn^j . Hence for a 

130 GeV gamma-ray line, the DM mass preferred is 
^142 GeV for annihilation to Z7 and a DM mass of 
^155 GeV is preferred for annihilation to H-f. Given 
the self annihilation cross section (av)~ ll presented in 
Ref. [14] , we convert them to (av) z / Hl by following 
the prescription given in Ref. [126]. For DM self- 
annihilation to Z7 or Hj, the relation between DM mass 
and the gamma-ray line is given by m x = (1/2) (1 + 

y/l + m%i h l 'E^)E 1 , and it follows from kinematic 

considerations that if the limits of (cxi>) 77 are given, the 
corresponding limits for (av) z/H-y is given by (av) z/H-y — 

(1/2) (l + y/l + m 2 z/h /E^j (crw) 77 . We take the upper 

and lower limits on (<jv) 77 for the 130 GeV gamma-ray 
line from the Region 4 of the SOURCE class events as 
presented in Ref. [14]. Using the other regions and the 
ULTRACLEAN class events gives similar limits and it 
will not change our conclusions. For the 408 MHz radio 
band, we only show the (av) z/H-y that is preferred by the 
130 GeV gamma-ray line for the NFW contracted profile. 

For both the DM self annihilation channels XX ~~ ^ %1 
and XX ~> Hj, we see that the contracted NFW profile 
gives the most constraining limit ((av) < 10 -27 cm 3 s _1 ), 



and in fact the sensitivity to the cross section is less 
than the total thermal relic cross section for both the 
self annihilation channels. The least constraining limit 
is obtained from the Einasto DM profile, as expected 
((av) < 10 -25 cm 3 s _1 ). If we assume that the modeling 
of the magnetic field near the GC black hole is correct, 
then this shows that the interpretation of the line signal 
near the GC for a contracted NFW profile is in mild 
tension with the radio data, provided the source of the 
gamma-ray line in the GC is due to the XX ~~ ^ ^1 an d 
XX ~^ H"f se lf annihilation channel. 



3. Sensitivity from a future measurement at 80 MHz 

The situation is expected to improve dramatically with 
future observation of the GC by LWA, LOFAR, and 
SKA. Although we present our future constraint from 
a radio flux measurement at 80 MHz near the GC, it 
is worth mentioning that strong constraints can also be 
obtained from measurement of the radio flux away from 
the GC. As mentioned earlier, ideally we expect the 
best measurement to come from a radio cold spot. The 
standard astrophysical background has to be modeled 
very carefully to reach the sensitivity as presented in this 
paper. 

To forecast the sensitivity, we very conservatively 
assume that LWA can reach a background subtracted flux 
density sensitivity 10 Jy at 80 MHz for a circular region 
of radius 2° around the GC [87]. The constraint in the 
(av)-m plane that can be derived from the radio flux 
measurement at 80 MHz for a circular region of radius 
2° is plotted the bottom panel of Fig. 7. We also show 
the (av) z/H-y preferred by the 130 GeV gamma-ray line 
by the green, red and blue shaded boxes for the NFW, 
Einasto and the contracted NFW DM profile respectively. 
We again use the Region 4 in the SOURCE class events 
and the prescription given in the previous section to draw 
these boxes. We did not draw these shaded boxes in the 
correct DM mass positions for clarity. 

Due to the superior flux sensitivity of LWA at these 
frequencies, we see that both the XX ^1 an d 
XX Hj self annihilation channel can be probed well 
below the total thermal relic cross section for all three 
considered DM profiles. In particular, for all the DM 
profiles considered, one can probe below the (av) ~ 
10~ 27 cm 3 s _1 cross sections required to explain the 
tentative 130 GeV signal. Thus, if the 130 GeV gamma- 
ray line turns out to be robust and originates from DM 
self annihilation, LWA has a good chance to search for 
the self annihilation channel giving rise to the line for 
the NFW, Einasto and the contracted NFW profile. Up 
to the uncertainty in the GC model, this remains, to 
our knowledge, the best probe for discerning the origin 
of the DM line independent of any particle physics DM 
model. Since LWA will reach this sensitivity over a large 
region of observation, the dependence of the constraint on 
the underlying DM profile is modest. We expect similar 
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limits can be obtained by the LOFAR collaboration 
as well. SKA is expected to further strengthen this 
constraint. 



V. SUMMARY AND OUTLOOK 

In this paper we have shown that existing radio 
data around the Galactic Center at 408 MHz marginally 
constrains the interpretation of the 130 GeV line in 
Fermi-LAT data in terms of DM self annihilation to 
Z^f or 7 with a cross section ~ 10 _27 cm 3 s _1 for a 
contracted NFW profile. For other frequencies or other 
DM density profiles the constraint is up to an order of 
magnitude weaker within the parameter ranges chosen 
by us. Future measurements made around the GC by 
LWA in the 80 MHz band can push the sensitivity to 
DM annihilating to gamma-ray lines down to (av) ~ 
10~ 28 cm 3 s -1 and enable a test of the above signal. 
Although the background needs to be known very well 
to achieve our quoted limits, these possibilities are, to 
the best of our knowledge, some of the most competitive 
ways to test for the nature of the DM that could have 
produced the tentative 130 GeV line signal. 

We have shown that these conclusions are fairly robust 
with respect to the assumptions on the magnetic field 
in the Galaxy, and the constraints do not weaken by 
more than an order of magnitude. The dependence 
on DM density profiles is somewhat more important, 
especially when the region of observation is small and 
closely centered on the GC. While the uncertainty in 



the astrophysical modeling of the GC does impact our 
results, we must emphasize that these constraints are 
completely model-independent from the particle physics 
perspective, because we have simply taken the electrons 
and positrons from the known decays of the Z or 
H produced in the DM annihilation to Z7 or Hj, 
respectively. 

We hope that these results will encourage radio 
astronomers, especially those at LWA, VLA-Low, 
LOFAR, and SKA, to observe the GC, model the 
astrophysical synchrotron backgrounds, and determine 
if there is any excess flux. Irrespective of whether the 
tentative 130 GeV gamma-ray line signal at Fermi-LAT 
is due to DM annihilation or not, this promises to deliver 
some of the strongest constraints on DM annihilation. 
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